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We perform an analysis of a binding energy difference called 5Vpn{N, Z) = —j \E{Z, N) — E{Z, N— 
2)~E{Z — 2, N)+E{Z—2, A^ — 2)] in the framework of a realistic nuclear model. It has been suggested 
that SVpn values provide a sensitive probe of nuclear structure, and it has been put forward as a 
primary motivation for the measurement of specific nuclear masses. Using the angular-momentum 
and particle-number projected generator coordinate method and the Skyrme interaction SLy4, we 
analyze the contribution brought to SVpn by static deformation and dynamic fluctuations around 
the mean-field ground state. Our method gives a good overall description of SVpn throughout the 
chart of nuclei with the exception of the anomaly related to the Wigner energy along the N = Z line. 
The main conclusions of our analysis of SVpn , which are at variance with its standard interpretation, 
are that (i) the structures seen in the systematics of SVp„ throughout the chart of nuclei can be 
easily explained combining a smooth background related to the symmetry energy and correlation 
energies due to deformation and collective fluctuations; (ii) the characteristic pattern of SVpn having 
a much larger size for nuclei that add only particles or only holes to a doubly-magic nucleus than 
for nuclei that add particles for one nucleon species and holes for the other is a trivial consequence 
of the asymmetric definition of SVpn, and not due to a the different structure of these nuclei; (iii) 
SVpn does not provide a very reliable indicator for structural changes; (iv) SVpn does not provide 
a reliable measure of the proton-neutron interaction in the nuclear EDF, neither of that between 
the last filled orbits, nor of the one summed over all orbits; (v) SVpn does not provide a conclusive 
benchmark for nuclear EDF methods that is superior or complementary to other mass filters such 
as two- nucleon separation energies or Q values. 

PACS numbers: 21.10.Dr; 21.30.Fc; 21.60.Jz 



I. INTRODUCTION 

Nuclear masses are measured today with an unprece- 
dented accuracy [H-Ill, in many cases better than a few 
keV. Such an accuracy is obtained not only for nuclei 
close to the stability line, but also for exotic ones with 
very short lifetimes. A recurrent question is how to take 
advantage of this major advance and how to use it to im- 
prove the theoretical description of nuclear ground states. 

The first possibility is to compare directly masses, or 
better binding energies, to theoretical predictions. Unfor- 
tunately, ab-initio methods based on a realistic nucleon- 
nucleon interaction are not available for systematic stud- 
ies of heavy nuclei. If they were, any disagreement with 
the experimental data would point to a deficiency of the 
interaction. When the many-body problem is not solved 
exactly (or, to be more precise, with a controlled numeri- 
cal accuracy), but with an effective model using effective 
degrees of freedom and an effective interaction, the link 
between data and nucleon-nucleon interaction is broken 
and a discrepancy between calculation and experiment 
can have its source in any ingredient of the model. 

The best available theoretical descriptions of masses 
[U, are not based on ab-initio methods. The three 
main models rely on very different ingredients. The mass 
formula of Duflo and Zuker Q does not make an explicit 
reference to a nucleon-nucleon interaction. Nevertheless, 



it assumes that there exist effective interactions smooth 
enough for Hartree-Fock calculations to be possible. The 
corresponding Hamiltonian is separated into monopole 
and multipolc terms that arc parameterized through scal- 
ing and symmetry arguments [ij- The macroscopic- 
microscopic approaches of MoUer et al. Q combine a 
finite-range liquid-drop or droplet model and shell ef- 
fects introduced through the Strutinsky shell correction 
method and a parameterized one-body potential. The 
main ingredient of the Hartree-Fock-Bogoliubov (HFB) 
mass formulae of Goricly et al. is an energy density func- 
tional (EDF) as widely used in self-consistent mean-field 
calculations. In a first variant, a Skyrme EDF is sup- 
plemented by empirical corrections for correlations that 
cannot be included in a mean field 01- In a second vari- 
ant, the same Gogny interaction is used to determine the 
mean field and quadrupole correlation effects beyond the 
mean field through a microscopic Bohr Hamiltonian 
The comparison of any of these models with data can 
hardly allow to extract general information about the 
nucleon-nucleon interaction. Comparison with results 
obtained using ab-initio methods and realistic interac- 
tions can be made through the idealized model of infinite 
nuclear matter. From such calculations, one can extract 
specific parameters, such as volume and symmetry energy 
coefficients for instance, corresponding to a liquid-drop 
formula (LDM) fitted to masses. Here also, the connec- 
tion between theory and experiment is ambiguous. The 
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liquid-drop model is justified by a Icptodcrmous expan- 
sion of the energy that cannot be expected to converge 
quickly even for the heaviest nuclei 0, [l^l ■ 

An apparently more appealing way to proceed is to re- 
late differences between binding energies of neighboring 
nuclei to specific features of nuclear models, in partic- 
ular to effective single-particle energies or effective two- 
body matrix elements . Models that fail to reproduce 
masses with a good accuracy are often more reliable for 
mass differences. The reason of this success is that mass 
residuals Mth — -/Wexpt for adjacent nuclei are not inde- 
pendent in a given mass model [l^ [l^ . 

This property has been used for a long time to asso- 
ciate one-nucleon separation energies with single-particle 
energies, or higher-order differences with pairing gaps. 
In particular, two-particle separation energies are impor- 
tant indicators of shell closures. One has to distinguish, 
however, between the use of mass filters as measures of 
specific model ingredients, and their use as signatures of 
structural changes. In particular, a mass filter cannot be 
expected to provide both simultaneously. One can test 
a model ingredient only when all nuclei entering a mass 
filter have the same structure, which becomes increas- 
ingly improbable with the number of nuclei involved. By 
contrast, the indication of a structural change (such as 
onset of deformation) by a mass filter often means that 
the fundamental assumptions made for its direct associ- 
ation with a feature of a model are violated. 

In a previous paper, we discussed the difficulties en- 
countered when trying to relate structures seen in the 
systematics of energy differences to features of the single- 
particle levels The changes in the large gaps ob- 
served in data for two-neutron or two-proton separation 
energies are often interpreted as due to the evolution of 
shell structure with TV and Z, and associated with the 
presence of strong residual tensor interactions [isj or a 
weakening of the spin-orbit interaction in neutron-rich 
nuclei [Hi [13 ■ We showed that experimental data can 
be explained in a coherent way within mean- field- based 
models as due to a combination of the slow modification 
of spherical single-particle spectra and the often rapid 
variation of collective correlation effects. 

In this paper, we perform a similar analysis for a mass 
filter that has become fashionable and has been identi- 
fied with the proton-neutron interaction in nuclei. Many 
recent experimental data have been used to interpret a 
difference between the (negative) binding energies of four 
even-even nuclei defined as 

SVpniN,Z) = -\[E{Z,N)-E{Z,N -2) 

-E{Z - 2, TV) + E{Z -2,N- 2)] 

(1) 

in terms of the effective interaction between the last oc- 
cupied neutron and proton orbits [l8l - l25| . It has been 
suggested that SVpn values provide a sensitive probe of 
nuclear structure, and it has been put forward as a pri- 
mary motivation for the measurement of specific nuclear 



masses [26|-|29[. A similar quantity has been analyzed in 
great detail by Janecke et al. [3^, but for nuclei differ- 
ing by one neutron and/or one proton only, which makes 
its interpretation more difficult. Indeed, the breaking 
of a pair in an nucleus with an odd number of particles 
modifies pairing correlations deeply and makes the struc- 
ture of its wave function significantly different from that 
of its even neighbors. By comparing only nuclei with an 
even number of neutrons and protons, one can hope that 
the assumption of a common mean field is better jus- 
tified. Differences between two consecutive SVpn values 
have also been proposed as a measure for the Wigner en- 
ergy that leads to an anomaly of binding energies around 
the N = Z line 

In Sect, ini we show how, in the framework of the 
Hartree-Fock (HF) method, SVpn can be related to the 
interaction between the last filled neutron and proton or- 
bits when making the same assumptions as those used to 
derive Koopman's theorem [s^ . We show that the intro- 
duction of pairing correlations and density dependencies 
complicate the relation even when making oversimplify- 
ing assumptions about the evolution of wave functions 
with N and Z. We discuss the main effects that make 
any direct identification of SVpn with a proton-neutron 
interaction doubtful. In Sect. IIVI we present results 
obtained from calculations using a realistic microscopic 
model using the Skyrmc energy density functional SLy4. 
We demonstrate how the successive inclusion of correla- 
tions from spherical to deformed mean-field calculations 
and further to symmetry restoration and configuration 
mixing permits to improve at each step the agreement 
with the experimental data, while, at the same time, los- 
ing in an ever increasing manner the simple interpreta- 
tion of SVpn- Section fVl summarizes our findings. 



II. ANALYSIS OF SVr,, 



Surprisingly, the abundant literature discussing the 
relevance of SVpn contains only very few analyses of 
its relation with the proton-neutron interaction in non- 
schematic models. Exceptions are the shell-model study 
of Heyde et al. 3J] , who underline that a simple interpre- 



tation of SVpn can only be given when there is just one 
dominant orbital for protons and neutrons each, and the 
nuclear DFT study by Stoitsov et al. [l^], who, however, 
focus their analysis on the overall excellent reproduction 
of data for SVpn with their model, rather than on the 
ingredients of the model that contribute to it. Below, 
we review the assumptions to be made to relate SVpn to 
the effective proton-neutron interaction in finite nuclei, 
and discuss their validity in the context of realistic nu- 
clear models based on the self-consistent mean field and 
taking the entire space of occupied single-particle orbits 
into account. We start with a simple two-body Hamilto- 
nian and discuss the corresponding energy in the context 
of HF (without pairing) and HFB (with pairing) . Then, 
we generalize the discussion to a more realistic effective 
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interaction that includes three-body forces or density de- 
pendencies. 



A. Frozen HF with two-body interaction 

Let us start from a Hamiltonian consisting of a kinetic 
energy term and an antisymmetrized two-body interac- 
tion: 



1 



E 



(2) 



When hmiting the A^-body wave function to a single 
Slater determinant, the minimum value of the e nerg y is 
obtained by solving Hartree-Fock (HF) equations [3J]- In 
this case the energy for a nucleus consisting of N neu- 
trons and Z protons is given by 
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where are the matrix elements of the kinetic en- 

— TV Z 

ergy operator and w„^/„„/ of the two-body interaction 
calculated in the single-particle basis that solves the HF 
equations. We have added superscripts N, Z to these ma- 
trix elements to recall that the HF equations are solved 
self-consistently and that, in general, the wave functions 
differ for each combination of iV and Z values. 

Let us assume for the moment that the HF single- 
particle basis is identical for the four nuclei with (A^, Z), 
(iV - 2, Z), {N, Z - 2), and (iV - 2, Z - 2). This "frozen 
HF" approximation leads to 



<5y,r'-— (iv,z) 
1 



— 4 \yN-l,Z-\,N-\,Z-\^ Vj^-x^z.N-X.Z 

'\-VN,Z~l.N.Z-l + VN,Z,N,z) ■ (4) 

The superscripts have been dropped, as the mean field is 
supposed to be the same for the four nuclei. ^ 

A further simplification can be obtained by noting 
that, for even N and Z, the two valence neutrons (in- 
dices N—1 and N) and protons (indices Z — 1 and Z) 



^ This assumption, however, cannot be expected to be valid for 
all nuclei. Each nucleus enters the expression for SVpn for four 
different nuclei, such that ultimately all nuclei had to have the 
same mean field for this assumption to be fulfilled. 



occupy time-reversed orbits in the HF solution, and that 
the matrix elements are equal two by two 



SV^/'"—{N,Z) = -^{vN,z-i,N,z-i 



VN,Z.N.Z) 



(5) 



Even in this simple case, the final result is not a single 
matrix element, but a combination of matrix elements 
between the valence particles. 

This derivation requires the same assumptions as those 
made to derive Koopman's theorem [33| which relates 
single-particle energies and one-nucleon separation ener- 
gies. Koopman's theorem, however, is known to have a 
very limited validity in nuclear physics, cf. the discussion 
in Ref. \\M and references therein. 



B. Frozen HF+BCS and HFB with two-body 
interaction 

It is obvious that the assumptions made in the previous 
section will rarely be justified, even approximately. Let 
us first examine the consequence the partial occupation of 
single-particle levels due to pairing correlations. We con- 
sider only pure proton-proton and neutron-neutron pair- 
ing, which can be justified for nuclei with N sufficiently 
different from Z . In the presence of proton- neutron pair- 
ing, the many-body state could not be written as the 
direct product of a proton and a neutron BCS state, and 
the energy would not be separable into proton and neu- 
tron components anymore. 

The HF+BCS or HFB expectation value of a two-body 
Hamiltonian ([2]) for a nucleus with N neutrons and Z 
protons evaluated in the canonical basis is given by 

n p 

I 1 \ " -N,Z 1 2 
+ '^nn'nn' '"n,N ""n' ,N 
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+ E 

n,n'>0 

+ E ^PPP'P' ^P^^ ^P^^ V,z V,z • (6) 

p,p'>Q 

The Ufc and Vk are real occupation amplitudes with 

Uk =UT.> 0, Vk = 



and u| 



1. Indices k and k 



refer to conjugate states, which for the ground states of 
even-even nuclei are connected by the time-reversal oper- 
ator. We use the usual convention where summation over 
all indices indicates summation over all k and k, whereas 
summation over "positive" indices means that the sums 
arc over all k, but not their conjugate levels k. The sec- 
ond index of the occupation amplitudes recalls that these 
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are occupation numbers for a nucleus with '^nN^^ 
neutrons and X]p^p,z = ^ protons. Except for the last 
two terms, the summation runs over positive and nega- 
tive values of n and p. 

To derive a simple expression for 6Vpn, a similar as- 
sumption as in the frozen HF case has to be made, namely 
that the canonical single-particle basis is the same for all 
four nuclei entering a given SVpn{N, Z). In this case, the 
matrix elements tkk and Vkk'kk' do not depend on N and 
Z . In addition, one has to assume that the solution of 
the HFB equations for neutron states docs not depend 
on the number of protons and conversely. This leads to 



rT^HFBf„,„: 

^ ^ pn 



where we introduced the shorthand Au^ 



n.N 



— -,2 



(7) 



_2 for the change of neutron occupation numbers 



when removing two neutrons, ^„ Ai 



2, and its 



homologue Avp ^ 



Vp z-2 for protons. 



C. Frozen HF with three-body or 
density-dependent interactions 

The forms derived in the previous sections were based 
on a two-body Hamiltonian. However, interactions de- 
rived from first principles contain at least three-body 
forces [1^, which are crucial to perform ab-initio calcu- 
lations with some predictive power. In a similar manner, 
the more phenomenological EDF-bascd methods must 
include terms of higher order than a two-body interac- 
tion, which is usually done by the inclusion in the EDF 
of terms with density dependencies. In particular, the 
saturation of nuclear matter cannot be satisfactorily de- 
scribed without taking into account these terms [36|. 

A three-body force adds a term 
J6^^,3,k,l,m,n^ijkl7nnala]ala,,a„^al to the Hamilto- 
nian In the HF approximation, its contribution to 
the binding energy is given by: 
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(8) 

p.p' ,p"—i 

using a notation analogous to that of Eq. In the 

frozen HF approximation, a Hamiltonian including two- 
body and three-body forces leads to 
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(9) 



The first three lines represent the two-body and three- 
body interactions between the last two neutrons and the 
last two protons. The last two lines of Eq. how- 
ever, contain a sum over all other nucleons, which is in- 
compatible with the interpretation of SVpn {N, Z) as the 
interaction between the last two neutrons and protons. 

Most energy density functional constructed for 
self-consistent mean-field calculations include a non- 
linear density-dependent two-body term of the form 

1 J2i,j.m,7i '^'ijki f{Pn + Pp) ala''jana,n, where p„ and pp are 
the local densities of neutrons and protons. Popular ex- 
amples for f{pn + Pp) range from the simple non-intcgcr 
powers of the total density f{x) — x", such as used with 
most Skyrmc and Gogny interactions, to the very elabo- 
rate density dependencies /(x) ^ [l-|-6(a;-|-(i)^]/[H-c(a;-|- 
d)^] used in modern density-dependent relativistic mean- 
field models [13, ■ Any density dependence that is not 
a simple polynomial of the total density gives rise to SVpn 
values that cannot be separated into neutron and proton 
contributions, even within the frozen HF approximation. 

We will not give explicit expressions for the HFB case 
with three-body forces or density-dependent terms. It 
should be obvious by now that even making the assump- 
tion of a frozen common canonical basis will lead to a 
lengthy and complicated expression for SVpn{N, Z) that 
does not allow for an intuitive interpretation. 



D. Discussion 

Even when making the drastic approximation that the 
four nuclei entering the calculation of 6Vpn{N, Z) can be 
described by the same mean field, the expression that is 
obtained for realistic models includes a summation over 
all single-particle levels. ^ Furthermore, the occupation 



^ ^ pn 



'{N,Z) 



A similar observation has been made by Van Isacker et al. [3ll | 
who quote that "&Vpn{N,Z) is an average neutron-proton in- 
teraction over the last fcvir nucleons" , but without giving any 
reference. 
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of the levels around the Fermi energy is affected by the 
addition or the removal of nucleons and the contribu- 
tion of each single-particle level to SV^^^'"^"" {N, Z) is 
weighted with the difference of its occupation between 
the nuclei. It is therefore doubtful that 6Vpn allows to 
isolate an empirical interaction between the last neutron 
and proton orbitals, as claimed in Ref. p^ . 

Moreover, one might wonder whether the approxi- 
mation of a frozen canonical basis necessary to derive 
Eqns. ([5]) and (O is ever satisfied. To change the num- 
ber of neutrons or protons by two induces rearrangement 
and polarization effects that modify the single-particle 
wave functions for both kinds of nucleons. Even if these 
effects are most often small, it should not be forgotten 
that SVpn {N, Z) is a tiny fraction of the total binding 
energy only, ranging from 10~^ in light nuclei to 10" 
in heavy ones. Therefore, even small rearrangement ef- 
fects can have a large impact on the values obtained for 
SVpn {N, Z) . We will analyze the validity of the frozen 
basis assumption for a few selected cases in Sect. IIVCI 
below. Finally, a self-consistent mean-field description of 
a nucleus provides a reasonable first approximation but it 
neglects correlations beyond the mean field that also con- 
tribute to the binding energy on the MeV scale [s^, li^l ■ 
These correlations cannot be cast in a simple form in- 
volving only the interaction between a few particles, and 
they also destroy the simple relation between SVpn and 
proton-neutron matrix elements. 

In the remaining part of this article, we investigate 
the importance of self-consistency, deformation, pairing, 
and configuration mixing for the description of data for 
SVpn We also analyze to which extend SVpn values can be 
identified with the effective proton-neutron interaction. 



III. THE MODELS 
A. The beyond-mean-field model 

Our method used to calculate binding energies for the 
ground states of even-even nuclei is described in detail 
in Refs. [H, H^. In our analysis, we use the energies 
as tabulated in [ilj, and we added a few nuclei in the 
vicinity of ^"^Pb. As effective interaction, we employ 
the SLy4 parameterization of the Skyrme energy density 
functional (42| for the mean-field channel in connection 
with a density-dependent zero-range pairing interaction. 

Starting from a set of mean-field calculations includ- 
ing a constraint on the axial quadrupole moment, two 
kinds of correlations beyond the mean-field are intro- 
duced. First, the deformed wave functions are projected 
on both fixed particle numbers and on angular momen- 
tum J = 0. In a collective model terminology, these 
correlations would be called rotational correlations. A 
second step of our method consists in the mixing of 
projected wave functions with different intrinsic axial 
quadrupole moment of the underlying mean-field state in 
a generator coordinate method (GCM). In the language 



of collective models, this corresponds to a vibrational cor- 
rection. The final wave function has the form 

\JMv) = ^ f.jAct)PimPNPz\q) ■ (10) 
q 

The ket \q) is a (paired) self-consistent mean-field state of 
axial quadrupole deformation q. The operators P/v, Pz 
and P^g project out the component with the particle 
numbers and angular momentum quantum number we 
are interested in. The weights /j,i/((?) defining the mixing 
of the projected wave functions with respect to q are 
obtained by variation of the total energy. 

We stress that there are no assumptions made in the 
model about the amplitude of the quadrupole fiuctua- 
tions introduced into the calculations. Depending on the 
structure of a nucleus, this amplitude either corresponds 
to a small vibration around a pronounced minimum, to 
a large-amplitude motion in a soft and wide potential 
well, or to the mixing of several states around coexisting 
minima in the deformation energy surface. 

In the following, we will compare results obtained from 
the energies determined using three wave functions that 
successively add quadrupole correlations: 

1. self-consistent spherical mean- field states = 0); 

2. the self-consistent mean-field minimum in the space 
of axial rcfiection-symmetric deformations |(7min), 
which might be spherical; 

3. the ground state obtained after configuration mix- 
ing of J = projected axial quadrupole. We refer 
to these wave functions in the following as projected 
GCM. The energy gained through these correla- 
tions will be called beyond-mean-field correlation 
energy in what follows. 

In each of these cases, the wave functions are projected 
on particle number. 



B. Liquid drop model 

The Strutinsky theorem [i^, |3] allows to decompose 
the binding energy into a "macroscopic" liquid drop part 
and a microscopic "shell correction" . In this picture, 
the macroscopic energy defined through the liquid drop 
model varies smoothly with A'' and without any cor- 
relation energies from deformation, shell effects, or fluc- 
tuations in collective degrees of freedom. It constitutes a 
reference with respect to which one can put into evidence 
all quantum effects. 

In some of the figures below we show macroscopic ener- 
gies calculated from a liquid-drop model whose parame- 
ters have been adjusted to reproduce the average binding 
energies of spherical nuclei calculated with the Skyrme 
interaction SLy4 (T^]. Besides the standard volume, vol- 
ume symmetry, surface and (direct) Coulomb terms, the 
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FIG. 1: (Color online) Maps of SVp„; From top to bottom: 
calculated with a spherical liquid drop formula having the 
average properties of the SLy4 interaction, by self-consistent 
calculations with SLy4 obtained assuming spherical nuclei, 
allowing for (axially) deformed shapes, and derived from J = 
projected configuration mixing calculations. The bottom 
panel shows experimental data. 



macroscopic model comprises surface symmetry, curva- 
ture and Coulomb exchange terms 

+(asurf + ^surf.sym ^ 
3 Z^e^ 3 



2/3 ^4/3g2 



where A = N + Z and / = 



N-Z 
N+Z ■ 



(11) 

The radius constant 



ro entering the Coulomb energies is determined from the 
nuclear matter saturation density po of SLy4 as Tq = 
3/(4^Po). 

The by far dominating contribution to 5Vpn comes 
from the volume and surface symmetry energies [l^l 



pn 



a, 



surf,sym 



(12) 



The global A ^ scaling factor in this expression originates 

from the denominator of the P A~ ^^^J factor in the 
symmetry and surface symmetry energy terms which do 
not cancel out in SVpn- 

There are two contributions to this term which have 
the same scaling [4^ : the first one is the difference in ki- 
netic energy between protons and neutrons that fill sep- 
arate potential wells, and the other the isovector part of 
the nucleon-nucleon interaction. The latter has a shorter 
range than the average distance between nucleons, such 
that in a semiclassical approximation it acts between 
nearest neighbors only, leading to the characteristic A 
and A^/"^ scaling of terms in Eq. p^ . The contribution 
of all other terms in Eq. ([TT]) to SVpn is not exactly zero, 
but it is too small to be resolved in the plots shown be- 
low. A standard liquid-drop model obtained with a "best 
fit" to experimental masses gives SVpn values that are 
systematically larger than those obtained from Eq. (|lip . 
mainly because the volume symmetry coefficient Osym has 
a slightly larger value than the one determined from the 
SLy4 interaction. 



IV. RESULTS 

A. Global behavior of SVpn 

The binding energies of Refs. [s^ and tabulated 
in [4l] | cover the region of even-even nuclei heavier than 
^^O for which experimental data are available, plus a few 
additional nuclei around doubly-magic systems. For the 
present study, we calculated a few extra nuclei around 
208p^^ Values obtained for SVpn with this sample of nu- 
clei are plotted as maps in Fig. [T] For a better resolution 
of the local fluctuations, the same data are plotted for 
isotopic chains as a function of the number of neutrons 
in Fig. [5J and for isotonic chains as a function of proton 
number in Fig. [31 
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The spherical macroscopic values are given each time 
in the top panel. All nuclei between the drip lines are 
represented for the LDM results in Fig. [U whereas in 
Figs. [2] and [3] results are restricted to the same set of 
nuclei shown in the other panels. The macroscopic SVpn 
values exhibit a regular smooth pattern and fall off with 
~ 1 /A. The slope of the decrease is related to the symme- 
try and surface symmetry energy coefficients of the EDF, 
Eq. p^ . This smooth systematic decrease of the SVpn 
values with increasing A has been sometimes interpreted 
as a result of "the gradual decrease in valence proton and 
neutron orbital overlaps due to the occupancy of shells 
of different average radii" [H, [22|. This is, at best, a 
mo del- dependent statement that cannot be translated to 
methods that calculate the energy from the interaction 
between all occupied particles. In particular, at no point 
in the derivation of the LDM expression p2)) has one to 
consider the form of the single-particle wave functions 
and their overlaps. Instead, it is only assumed that all 
occupied single-particle wave functions add up to the sat- 
uration density inside the nucleus. 

The smooth trend of the macroscopic calculation is 
still apparent in the spherical self-consistent mean-field 
results. Some deviations appear, however, which are re- 
lated to the magic numbers at 20, 28, 50, 82 and 126. 
For nuclides just below these shell closures, the spherical 
mean-field results are slightly larger than the LDM ones, 
whereas they take slightly smaller values for nuclei just 
above. As a consequence, SVpn values do not fall off con- 
tinuously with A, but form sheets separated by the shell 
closures. 

Relaxing the constraint of spherical symmetry strongly 
modifies the behavior of dVpn by giving rise to rapid fluc- 
tuations around the smooth trend, with an amplitude of 
up to 200 keV. This change can be directly related to the 
effect of deformations on binding energies. The variation 
of quadrupole deformation and of the associated energy 
gain with N and Z over the entire nuclear chart have 
been presented in Figs. 9 and 16 of Ref. The energy 
gain due to deformation can reach more than 20 MeV 
and can vary rapidly from one nucleus to the other. Any 
mismatch in the evolution of deformation energy between 
the four nuclei entering SVpn can dramatically change its 
value. 

To illustrate the impact of deformation and correla- 
tions on SVpn, the quadrupole deformation /32, deforma- 
tion and correlation energies for nuclei heavier than ^'^^Sn 
are presented in Fig.m As expected, the absolute value of 
the deformation energy increases flrst slowly when mov- 
ing away from the proton and neutron shell closures and 
then more rapidly until it peaks at almost 18 McV in 
the rare-earth region and above 20 MeV for actinides. 
As indicated in the inset of the flgure, SVpn is defined 
as the sum of the energies of the nuclei on the diagonal 
minus the sum of the energies of the nuclei on the anti- 
diagonal. Figure U gives an intuitive illustration of how 
the deformation and beyond-mcan-field correlation ener- 
gies contribute to SVpn- Along a line going from ^■^^Sn to 
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FIG. 2: (Color online) Same data as Fig. \T\ but plotted for 
isotopic chains as a function of neutron number. 



^"^^Pb and beyond, the deformation energy varies rapidly 
and nonlincarly and brings a very large contribution to 
SVpn {N, Z) . On the contrary, for nuclei located close to 
a spherical shell closure for one nucleon species and mid- 
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FIG. 3; (Color online) Same data as Fig. [1] but plotted for 
isotonic chains as a function of proton number. 



shell for the other, the lines of equal deformation and 
correlation energy are nearly parallel to the N ov Z axis, 
leading to a pairwise cancelation of similar deformation 
energies. The fine structure of this cancelation depends 
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FIG. 4: Effect of deformation and correlation energies on 
SVpn values for heavy nuclei. Top to bottom: (i) map of the 
dimensionless quadrupole deformation /32 of the mean-field 
ground state, (ii) contour plot of the static deformation en- 
ergy, (iii) contour plot of the sum of the static deformation 
and dynamical beyond-mean-field correlation energies, and 
(iv) map of the contribution of deformation and dynamical 
correlation energy to 5Vpn- The total 5Vpn value is then ob- 
tained adding the spherical mean-field value. The inset gives 
a reminder of the relative signs of the four contributions to 
SVpn, their distance being drawn on the same scale as the one 
used in the contour plot. 




contribution of deformation and 
correlation energies to dV.,,; 



■ — I 



0.1 
0.05 
z\ 0.0 (MeV) 

1-0.05 
-0.1 



9 



of course on the deformed shell structure of the four nu- 
clei entering a given dVpn {N, Z) value and leads to an 
erratic behavior of 5Vpn, particularly visible when drawn 
for isotopic or isotonic chains as in Figs. [5] and [31 The 
same behavior is also seen in the experimental data. It 
is worthwhile to mention that there is no direct relation 
between the size of the deformation and the deformation 
energy, nor between the sign of the deformation, prolate 
or oblate, and the deformation energy. In particular, the 
transition between oblate and prolate shapes in a region 
of shape coexistence around the neutron deficient Pb iso- 
topes does not leave obvious traces in the ground-state 
deformation energy, hence, in the calculated SVpn val- 
ues. As one may expect, introducing beyond-mean-field 
correlations evens out the effect of static deformation. 
The effect of configuration mixing is, indeed, a spread- 
ing of the ground-state wave function around the mean- 
field minimum and a mixing of coexisting shapes. The 
bcyond-mean-ficld correlation energy varies rapidly only 
around shell closures and has its largest impact on the 
5Vpn{N,Z) values in these regions. We will analyze its 
impact in more detail for selected nuclei below. 

The very rapidly varying behavior of 5Vpn around the 
N = Z line that sticks out in the experimental data for 
light nuclei in Figs. [TJ [2] and [3] is not reproduced by any 
of our calculations. This anomaly is due to the Wigner 
energy [sil IH, IH, \^ , whose origin is not described by 
present-day EDF models, see Refs. [H, HI] for further 
discussion of this deficiency. 

In the literature one cannot find, however, a unique 
definition of the Wigner energy. Sometimes this notion 
is used for an anomalous additional contribution to the 
binding energy of the T = N = Z member of an iso- 
baric multiplet compared to the {N — Z)^ extrapolation 
from the other isobars H^, but more often the Wigner 
energy denotes a contribution to the binding energy that 
is linear in jT^j = — Z|. Such a term arises, for exam- 
ple, from a Hamiltonian that is invariant under Wigner's 
SU{4) symmetry, E ~ T{T + 4) = E,y„, + Ewisnc- 
When using this second concept of a Wigner energy, 
the anomaly of binding energies at iV = Z is the con- 
sequence of the Wigner energy -Ewigncr having a discon- 
tinuity in its derivative at = 0, and not of an ad- 
ditional binding of N = Z nuclei. There are many rea- 
sons whyWigner's SU{A) symmetry is not realized in nu- 
clei [sTI. Is^ lia. lirl. [sof ■ which suppresses the linear term 
in \N — Z\ compared to the quadratic one. Still, traces of 
such a linear term are implicitly contained in all realis- 
tic shell- model calculations [s^, [3 El, and explicitly 
in the Duflo-Zukcr mass formula [5|- It is noteworthy 
that the macroscopic-microscopic mass models and 
the Skyrme-HFB mass formulae also contain explicit 
phenomenological corrections for the Wigner energy that, 
in fact, combine both of the concepts of Wigner energy 
mentioned above, see also the discussion in Ref. Sev- 
eral mass differences have been put forward as indicators 
or even measures of the Wigner energy, most prominently 
double- /3-decay Q values [iy, 113] or a difference between 
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FIG. 5: (Color online) Top to bottom: Intrinsic deforma- 
tion of the mean-field ground state, 1/4 of the deformation 
energy, 1 /4 of the beyond- mean-field correlation energy, and 
SVpn values obtained from spherical mean-field calculations, 
deformed mean-field calculations, J = projected GCM cal- 
culations and experimental data for the chain of Sn {Z = 50) 
isotopes. For the deformation, deformation energy and cor- 
relation energy the values for Cd {Z = 48) isotopes are also 
shown. 



three different SVon values 



It has to be stressed. 



however, that SVpn itself is not a measure of the Wigner 
energy. 



B. Selected chains of nuclei 

1. Isotopic chains of magic nuclei 

To demonstrate how 6Vpn is build up from different 
types of correlations, let us examine now its evolution 
along cuts through the N~Z plane. We first look at the 
two isotopic chains of Sn and Pb, corresponding to closed 
proton shells, in Figs. [5] and [6] We showed in Ref. 
for the Z = 50 Sn isotopic chain that static deformation 
and dynamical correlation energies are key ingredients 
in reproducing the two-proton separation energies across 
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FIG. 6: (Color online) The same as Fig.O but for the chain of 
Pb {Z = 82) isotopes. The oblate ground-state deformation 
found for Hg isotopes is confirmed by experiment below = 
122. 



the Z ~ 50 shell and in explaining in an intuitive way 
the mutually enhanced magicity around ^^^Sn. 

The results for Sn isotopes are displayed in Fig. [S] and 
those for Pb in Fig.|6l The dimensionless deformation pa- 
rameter ^2^^'^ is related to the intrinsic mass quadrupole 
moment of the self-consistent mean-field wave functions 
(q|(2z^ - .T^ - y'^)\q) as 
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quantity entering 6Vpn, Eq. ([T|) is plotted on the figure. 
The magnitude of the contributions of deformation and 
correlations to 5Vpn{Z, N) is directly related to the dif- 
ference in slopes of the curves for Z or Z — 2 for a given 
N. The largest contributions are obtained when one of 
the slopes is much steeper than the other. 

The ground-state configurations obtained for most Sn 
and Pb isotopes are spherical; some mid-shell isotopes are 
slightly deformed. However, the energy gain due to de- 
formation in those cases is small, smaller than 200 keV, 
and it originates from a deformed minimum nearly de- 
generate with the spherical configuration. By contrast, 
the deformation and gain in deformation energy for the 
ground states of the non-magic Z — 2 isotopic chains can 
be large and vary rapidly for some neutron numbers. For 
those cases, the contribution of the deformation energy 
to 6Vpn is large. This clearly indicates that one cannot 
assume to describe all the four nuclei entering SVpn by a 
common mean-field, even for closed-shell nuclei. 

The correlation energy is larger for all Cd and Hg iso- 
topes than for Sn and Pb nuclei with same N. However, 
the slopes of the Z and Z — 2 curves differ significantly 
for a few isotopes only, and its contribution to SVpn is 
large only for these neutron numbers. Nevertheless, the 
beyond-mean-field correlations level out the rapidly fluc- 
tuating effect of static deformations. Their contribution 
to the binding energy plays also a key role in the de- 
scription of the two-proton se par ation energies across the 
Z = 50 and Z = 82 shell [li,|4|. 

The agreement between the experimental data and the 
results of the beyond-mean-field calculation is very sat- 
isfactory for the Sn and Pb isotopic chains, as can be 
seen in the bottom panel of the Figs. [S] and [51 In par- 
ticular, only the latter calculation is able to describe the 
rise of SVpn up to iV = 126 and its sudden drop beyond. 
The rapid variation of SVpn around ^'^^Sn and ^osp^^ jg 
mainly due to the onset of substantial beyond-mean-field 
correlations around doubly-magic nuclei. This scenario 
is much more involved than the proton-neutron interac- 
tion between the valence orbitals invoked in Ref. [2^. A 
detailed analysis of the contributions to the SVpn value 
of 2"8pb is given in Sect. HVCl 
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with Rq = 1.2 A^/'^ fni. In the three upper panels of the 
figures, the results for isotopic chains with two protons 
less are also displayed, such that the values for all four 
nuclei entering SVpn are given in the same plot. The 
contribution of deformation and correlation energy to a 
given SVpn {N, Z) value can be extracted from the plot by 
first taking the difference between the energies for N and 
A'' — 2 on the curves for Z and Z — 2, and then subtract- 
ing the value for Z — 2 from the one for Z. To facilitate 
this comparison, the deformation and correlation ener- 
gies have been multiplied by a factor 1/4, so that the 



2. Onset of deformation in rare-earth nuclei 

Let us now analyze the isotopic chains of Ba {Z = 56), 
Nd (Z = 60), and Gd {Z = 64), which cover a region 
of nuclei with a large variation of deformation on both 
sides of the spherical = 82 shell closure. Results of our 
calculations are compared with the experimental data in 
Figs. [71 m and [HI Intrinsic deformations calculated with 
Eq. are compared to values taken from Refs. [sil . 
|52| that are determined from experimental B{E2) values 
making the assumption of a rigid axial rotor 
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FIG. 7: (Color online) Top to bottom: Intrinsic deformation 
of the mean- field eround state compared with experimental 
data taken from [52|, 1/4 of the deformation energy, 1/4 of 
the beyond-mean-field correlation energy, and 5Vpn values ob- 
tained from the macroscopic model, spherical mean-field cal- 
culations, deformed mean-field calculations, J = projected 
GCM calculations and experimental data for the chain of Ba 
{Z = 56) isotopes. For the deformation, deformation energy 
and correlation energy the values for Xe {Z = 54) isotopes 
are shown as well (see text). 



with i?o = 1.2A^/'^fm. For well-deformed nuclei, theo- 
retical ([T3| and experimental (fT4|) values are in excellent 
agreement. Around spherical shell closures, however, the 
lowest 2+ state is dominated either by non-collective two- 
quasiparticle configurations or by fluctuations in collec- 
tive degrees of freedom, neither of which can be described 
by the mean-field ground state. 

For none of these three isotopic chains, the spherical 
mean-field result for SVpn does show any structure except 
for a tiny drop at the N = 82 shell closure that becomes 
rapidly smaller with increasing proton number. Besides 
that, the spherical mean- field values remain very close to 
the macroscopic ones for all three chains of nuclei. The 



FIG. 8: (Color online) The same as Fig. [71 but for the chain 
of Nd {Z = 60) isotopes. 



only isotopes to remain spherical when deformations are 
allowed, arc those with TV = 82. All heavier isotopes are 
prolate, with a very similar variation of deformation as a 
function of N for all of the three isotopic chains. Lighter 
isotopes are prolate for Ba and Nd, and oblate for Gd. Al- 
though the deformation varies with A'^ in a rather similar 
way for the three chains, the effects of deformation and 
beyond-mean-field correlations on SVpn are different. De- 
formation and correlations do not bring very large contri- 
butions, but they induce a significant change of behavior 
of 6Vpn for Ba and Nd. Note also that the changes with 
respect to the spherical case bring theory closer to exper- 
iment with a very few exceptions. The contribution from 
deformation energy overcorrects the spherical result for 
SVpn , in particular by making SVpn smaller below A'^ = 82 
and larger above this value. The beyond-mean-field cor- 
relations straighten the curve and bring it very close to 
the data. 

Comparing the three isotopic chains, the largest devia- 
tion between the experimental and the macroscopic SVpn 
values is observed for the Ba {Z = 56) isotopes. For the 
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FIG. 9: (Color online) The same as Fig. [71 but for the chain 
of Gd [Z = 64) isotopes. 

isotopes below iV = 82 (with the possible exception of 
the lightest one N = 66), the experimental values are 
smaller, whereas above = 82 they are larger. The 
same overall behavior is also found for Nd {Z = 60), but 
with a smaller deviation from the macroscopic results. 
For Gd {Z = 64), the experimental data lie almost on a 
straight line, very close to the macroscopic results. 

The EDF models provide a simple explanation of these 
different behaviors. The three chains present a similar 
evolution as a function of N, going from deformed to 
spherical to deformed shapes again. However, looking to 
Fig. 01 one can see that the chains are located differently 
with respect to the center of the deformed region. The 
Ba {Z ~ 56) isotopes are situated at the lower end, where 
the deformation energy grows with at a very different 
rate for Z and Z — 2 and brings a large contribution to 
SVpn- The Gd {Z = 64) chain is close to the center of the 
deformed region where the deformation energy of adja- 
cent isotones grows synchronously. The SVpn values are 
not only unaffected by deformation in the Gd isotopes, it 
is also remarkable that the shape transition from a pro- 
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FIG. 10: (Color online) The same as Fig. [71 but for the chain 
of constant N — Z = 32. The nuclei for which data exist are 
expected to be prolate, whereas the systematics of rotational 
bands and radii in the Hg isotopic chain suggests that '^^''Hg is 
oblate in its ground state, in agreement with the calculation. 

late shape for Z = 62 to an oblate one for Z = 64 at 
N = 78 does not visibly affect the SVpn value obtained 
from the deformed mean-field calculation. These exam- 
ples indicate that SVpn cannot always be expected to be 
a sensitive indicator for changes in deformation. 

3. SVpn along lines of constant N + Z or N ~ Z 

Let us now examine the evolution of SVpn along other 
cuts through the chart of nuclei. In Fig. (TUl results for 
nuclei on a line of constant T, = N — Z ^ 32 are provided 
as a function oi A = N + Z. Available data start at the 
doubly-magic nucleus ^'^^Sn, cover the entire rare-earth 
region, and extend beyond Z = 82. Most of the nuclei 
along this chain are deformed, except for Z = 50 and 
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FIG. 11: (Color online) The same as Fig. □ but for the 
isobaric chains with A = N + Z = 132 and A — 168. 



82. The deformation energy takes its largest absolute 
value of around 16.8 MeV for ^^®Er at Z = 68. However, 
it is not the magnitude of the deformation energy that 
governs the size of its contribution to (5V^„, but how the 
curve for Tz = 32 diverges from those for Tz = 30 and 
34. Indeed, the contribution of the deformation energy to 
SVpn is obtained as the sum of the energies of two succes- 
sive points on the Tz = 32 line from which one subtracts 
the sum of the values on the = 30 and — 34 lines 
corresponding to the intermediate A value. The contribu- 
tions of deformation and beyond-mean-field correlations 
reinforce themselves for the lightest nuclei, leading to a 
rather irregular pattern, significantly different from the 
structureless spherical results. It is remarkable that (5V^„ 
values obtained from the beyond-mean-ficld calculation 
follow very closely the many irregularities of the experi- 
mental data. 



ing to fixed values of A, which are perpendicular to the 
TV - Z = 32 chain, are plotted in Fig. [iTl The first 
one corresponds to A = 132 and extends from the very 
deformed neutron-deficient ^'^^Nd to the doubly magic 
^■^^Sn. As in the case of the Ba and Nd isotopic chains, 
the different onset of deformation and correlation energy 
for Z and Z — 2 lowers (5V^„ relative to the macroscopic 
values below the = 82 shell closure, whereas it is en- 
hanced for ^■^^Sn. Note for this nucleus the significant 
differences between the deformed and the beyond-mean- 
field calculations. The change of behavior for A^ — Z = 32 
of SVpn with respect to the spherical and macroscopic 
values is very nicely described by the beyond-mean-ficld 
calculation. 

The second isobaric chain in Fig. [11] A = 168, almost 
follows the diagonal in Fig. S) All isobars are deformed 
and the deformation increases gradually when going from 
the very neutron-deficient ^^^Pt to ^^^Er, a nucleus lo- 
cated in the center of the deformed rare-earth region. 
The deformed mean-field and beyond-mean-field calcula- 
tions give very similar SV^n values and agree well with the 
data. For the lighter isobars, the (5V^„ values are smaller 
than the macroscopic ones, whereas for ^^^Er they sud- 
denly increase to values above. For the chains that we 
discussed up to now, the sudden increase of SV^n from be- 
low to above macroscopic values took place when crossing 
a spherical shell closure, i.e. with decreasing deformation. 
In the case of the A = 168 chain, the sudden increase of 
bVpn has its origin in the saturation of deformation en- 
ergy with increasing asymmetry. 

Again, the deviation of SV^n from the macroscopic 
value depends on the difference of increase in deformation 
energy in adjacent nuclei. 



^. Doubly-magic nuclei and mutually enhanced magicity 

In Figs.[5l[ini andim the SVpn value of ^^'^Sn sticks out 
as being larger than that of all surrounding nuclei. The 
same result is obtained for •^°*Pb. In both cases, these 
SVpn values are also much larger than the macroscopic 
trend. A similar singular behavior for doubly-magic nu- 
clei is also found with other mass filters, such as two- 
particle separation energies or Qa values: the value ob- 
tained for a doubly-magic nucleus is much larger than 
those of adjacent nuclei, including the semi-magic ones. 
This gives the impression that the shell closure of one 
nucleon species reinforces the magicity of the other, an 
effect sometimes called "mutually enhanced magicity" in 
the literature P, IH, HJ]. This effect is not described 
by pure mean-field models for which two-nuclcon sepa- 
ration energies or Qa values across a shell closure usu- 
ally show very little variation with the number of parti- 
cles of the other species, in stark contrast to the data. 
In Refs. [ij, HO, Hall it was shown that these filters are 
much better described when beyond-mean-field correla- 
tions are taken into account. The same result is found 



Results obtained for two nuclear chains correspond- here for SVpn ■ The beyond- mean-field correlation energy 
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is much smaller in a doubly-magic nucleus than in its 
neighbors. Its rapid variation gives a contribution to 
SVpn that pushes it to very large values in doubly-magic 
nuclei, up to twice as large as the average trend. 

At the same time, the pattern of the SVpn values 
changes for doubly-magic nuclei. For nuclei located ei- 
ther below or above the shell closures for both nuclcon 
species, the experimental SVpn tends to be larger than 
the average trend. In contrast, for nuclei where one nu- 
cleon species is below and the other above the respective 
shell closure, the experimental SVpn value tends to be 
smaller than the average trend. This behavior is very 
well illustrated in Figs. [3 [8] and [9] for the shell closures 
at Z = 50 and N = 82. 

In the literature, qualitative explanations have been 
proposed for this effect, based on the nature of the or- 
bitals filled by neutrons and protons. If the energies of 
both orbitals are larger or smaller {p — p and h — h cases) 
than that of the Fermi level, they are supposed to have 
a large overlap. On the contrary, if one of the energies is 
larger and the other smaller, [p ~ h and h — p case), this 
overlap is supposed to be small. The behavior of SVpn 
is then attributed to the differences between the over- 
laps ^ij. EDF calculations offer a different and more 
straightforward explanation. This effect results from the 
combination of a smooth macroscopic background and 
the contributions from deformation and beyond-mcan- 
field correlation energies. As can be seen in Fig. |4l their 
combined absolute value increases in all directions around 
a doubly-magic nucleus. Moreover, looking for instance 
at Figs. [5] and [HI one sees that this increase is nonlinear. 
The pattern that is observed for SVpn around doubly- 
magic nuclei is then a trivial consequence of the asym- 
metry of the relative signs of the four energies entering 
its definition. For example, for N = 118, the Z = 7A nu- 
cleus is located in such a way that the nonlinear increase 
of these contributions is pointing toward N — 2, Z ~ 2. 
For the same value of N but Z = 7% the iso-energy lines 
are nearly parallel to the A'^-axis and the contributions 
for a given Z-value nearly cancel out. 

Looking once more at Fig. 21 the picture on how SVpn 
is build up emerges clearly. Let us divide the map into 
rectangles delimited by the proton and neutron magic 
numbers. In any of these rectangles, deformation and 
correlation energies grow nonlinearly from small values 
along all borders to large ones in the middle. The result- 
ing map of deformation and correlation energies is highly 
symmetric and centered around the middle of the region. 
The definition of SVpn, Eq. ([Ij, however, is asymmetric. 
It is designed to probe the increase of the energy when 
going from the lower left to the upper right in the nu- 
clear chart under the assumption that it is superimposed 
on a background of like-particle interactions independent 
on the number of the other particle. The combined de- 
formation and dynamical correlation energy does rarely 
follow this anticipated pattern. As a consequence, one 
always obtains positive contributions to SVpn around the 
so-called p — p and h — h corners of the rectangle and 



negative values in the p — h and h — p corners. Close 
to the center of such a major-shell region, this trend is 
inverted when the deformation and correlation energies 
reach their maximum. This explanation of the pattern of 
SVpn around shell closures does not invoke any knowledge 
about the spatial structure of the single-particle orbits 
and their overlaps, and indicates also that the observed 
pattern of SVpn docs not necessarily signal stronger or 
weaker proton-neutron interactions in the four corners of 
a region of the nuclear chart between major shells. 

The observation that the appearance of enhanced SVpn 
values in the rare-earth region when going from ^"^^Sn to 
^°®Pb and beyond, is correlated to the line of A'vaicncc ~ 
Zvaionco has led the authors of Ref. [2^ to the speculation 
that this phenomenon might be due to a "mini-Wigner 
energy" of origin similar to the Wigner energy that leads 
to enhanced SVpn values along the N ~ Z nuclei. Our 
analysis makes this scenario very unlikely and offers a 
simpler explanation. First, we underline that our model 
does not give any trace of the Wigner energy and its con- 
tribution to SVpn at the N = Z line, as is the case for all 
present-day self-consistent mean- field models [H, [13, El] , 
meaning that the relevant physics is not contained in it. 
In contrast, our model does reproduce very well the en- 
hanced SVpn values along the A^vaience ~ ^valence huc in 
the rare-earth region. As explained above, their enhance- 
ment is a consequence of the onset of deformation and 
beyond-mean-ficld correlations when going away from a 
doubly-magic nucleus, which gives a positive contribution 
to SVpn in some direction and negative in other directions 
due to the asymmetric definition of SVpn ■ 



C. Detailed analysis of 5Vp„ for selected nuclei 

1. General comments 

The discussion above demonstrates that the rapid 
variation of the deformation energy and the beyond- 
mean-field correlation energy from symmetry restora- 
tion and shape mixing often gives large contributions to 
SVpn{N, Z). For nuclei away from the N ~ Z line, this 
variation is at the origin of almost all structures seen in 
the data. This also indicates that the assumption of a 
common single-particle basis made in Sect. |II]to obtain 
a simple expression for SVjjn in terms of proton-neutron 
matrix elements is rarely justified. When the structure of 
the four nuclei entering Eq. ([T]) is different, the question 
arises whether there are other terms in the energy func- 
tional than the proton-neutron interaction energy that 
contribute to SVpn- 

We have selected three representative nuclei for which 
we will decompose SVpn{N,Z) into contributions from 
the proton-proton, neutron-neutron and proton-neutron 
terms in the EDF. 
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TABLE I: Decomposition of SVpn into contributions com- 
ing from the different terms of the energy density func- 
tional for ^"^Pb, namely kinetic energy of neutrons and pro- 
tons, the neutron-neutron, proton-proton and proton-neutron 
parts of the Skyrme EDF, the neutron-neutron and proton- 
proton parts of the pairing functional and the proton-proton 
Coulomb EDF. We also give the sum of all terms and the 
experimental value. All energies are in MeV. 
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frozen HP 


HP 


HP-HBCS-HLN J 


= GCM 


kinetic n 


0.000 


-0.056 


-0.010 


-0.105 


kinetic p 


0.000 


0.012 


-0.044 


-0.132 


Skyrme nn 


0.025 


0.076 


-0.005 


0.160 


Skyrme pp 


0.017 


-0.008 


0.025 


0.275 


Skyrme pn 
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0.211 


0.218 


0.382 


pairing nn 






0.027 


-0.057 


pairing pp 






-0.005 


-0.081 


Coulomb 


0.000 


0.010 


0.020 


0.012 


total 


0.204 


0.245 


0.225 


0.457 



cycxp 



0.427 



2. 



*Pb 



The first nucleus ^"^Pb has been chosen for two rea- 
sons. First, we have seen in Fig. IH] that the contribution 
of the beyond-mean-field correlation energy to the 5Vpn 
value of this doubly-magic nucleus is particularly large. 
Second, the spherical ^osp^ presents a very favorable sit- 
uation to numerically test the frozen HF approximation, 
where the same set of single-particle wave functions is 
used to compute the energy of all four nuclei involved 
in the computation of SVpn- As discussed in Sect. Ill A| 
this approximation has to be made to establish the direct 
relation between 5Vpn and the two-body proton-neutron 
interaction. In fact, ^ospj^ jg Qj-^g ^j^g very few spherical 
nuclei for which such calculations can be performed. It 
requires that four neighboring nuclei have a closed shell 
configuration, which is possible only for N and Z values 
for which the orbitals below the Fermi level are Pi/2- 
or S1/2+ levels for both protons and neutrons. These 
conditions are met for ^°^Pb, with a v S1/2+ level below 
N = 126 and a tt pi/2- level below Z = 82. 

The results are presented in Table HI In the first 
column, the four nuclei entering 5Vpn have been cal- 
culated with the single-particle basis of ^"^Pb ("frozen 
HF" approximation), without readjustment of the basis 
for each nucleus and without pairing correlations. Self- 
consistency for ^"^Pb, ^"^Hg and ^""'Hg has been consid- 
ered for the results given in the second column ("HF"), 
and self consistency and pairing correlations treated with 
the BCS-f LN prescription have been taken into account 
for the values of the third column. Finally, the fourth 
column corresponds to the J = projected GCM calcu- 
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FIG. 12: (Color online) Difference between the values of the 
terms of the energy functional listed in Table [l] for ^"^Pb, 
2°'^Hg and ^'^^Hg and their value for 2°*Pb. The same four 
sets of calculations as in Table U are considered. All energies 
are in MeV. 



lation. 

We decompose the energy density functional into the 
kinetic energies of neutrons and protons (including the 
centre-of-mass correction), the neutron-neutron, proton- 
proton and proton-neutron parts of the Skyrme EDF that 
models the particle-hole part of the effective strong inter- 
action, the neutron-neutron and proton-proton parts of 
the pairing functional and the proton-proton Coulomb 
energy. The Skyrme and pairing functionals contain 
density-dependent terms. We interpret them as a 
density-dependence of the respective neutron-neutron, 
proton-proton and proton-neutron terms. This choice of 
decomposition is not unique, however. For further details 
about the functional, we refer to Refs. [36l[56j. 

In the frozen HF calculation, the sole contribution 
to SVpn comes from the Skyrme EDF. The neutron- 
proton terms give the largest contribution, although the 
neutron-neutron and proton-proton terms contribute by 
about 20% through their density-dependence. As soon as 
self-consistent wave functions arc used, the one-body con- 
tribution from the kinetic energy becomes large. This is 
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not surprising, as the kinetic energy provides a large con- 
tribution to the symmetry energy coefficient agym of the 
EDF [HB, [131 , which in turn dominates the global trend of 
SVpn, Eq. (fT^ . All other terms in the functional are mod- 
ified and can bring sizable contributions to 6Vpn ■ Pairing 
correlations change all contributions even further. The fi- 
nal value of SVpn from a self-consistent calculation results 
from a partial cancelation of several terms. The proton- 
neutron terms in the Skyrme functional are of the right 
order of magnitude and approach the final value of 6Vpn 
at a given level of approximation, but even for ^°*Pb, 
which is probably one of the most favorable cases for the 
frozen approximation, one can hardly conclude that SVpn 
is a valuable measure of proton- neutron interactions. 

More importantly, the spherical mean-field values are 
far from the experimental data and the correlations 
brought by symmetry restoration and configuration mix- 
ing arc large and crucial to obtain the correct value, as 
they increase SVpn by 0.232 MeV to almost twice the 
mean-field value. The decomposition of the J = pro- 
jected GCM results is listed in the fourth column of Ta- 
ble HI All terms become large in absolute value, but par- 
tially cancel each other. Their sum gives a value for SVpn 
close to the data. This is not entirely surprising as cor- 
relation energy is always gained from the compensation 
between a loss in kinetic energy and a gain in interaction 
energy. Table H] indicates that there are large contribu- 
tions to SVpn from proton-proton, neutron-neutron and 
neutron-proton terms and that none of them is dominant. 

Fig. [12] illustrates how the SVpn is build up from can- 
celations between the contributions of the four nuclei en- 
tering its definition. It shows the differences between 
the EDF terms defined in Table U for ^ospb, 206jjg and 
^"''^Hg and their value in ^°®Pb. The contribution of a 
given term to SVpn{^°^Ph) is the sum of the values for 
206pj-, (j-gj^ and ^^^Hg (green) minus the value for ^°*Hg 
(blue), divided by four. While in the frozen HF approx- 
imation, most terms cancel out, exactly self-consistency, 
pairing and collective quadrupole correlations lead to a 
much more complex situation. One can conclude from 
this analysis that each contribution to SVpn itself is the 
result of a partial compensation between changes in all 
four nuclei entering its definition. 



TABLE II: Decomposition of SVpr^ values for ^"^^Er and ^^^Nd 
into contributions coming from the different terms in the 
energy density functional in deformed self-consistent mean- 
field calculations, namely kinetic energy, the neutron-neutron, 
proton-proton and proton-neutron parts of the Skyrme EDF, 
pairing energy and the proton-proton Coulomb EDF. We also 
give the total contribution from beyond-mean-field correla- 
tions, the sum of all these terms and the experimental value. 
All energies are in MeV. 



term 


168 Er 


152Nd 


kinetic 
Skyrme nn 
Skyrme pp 
Skyrme np 
Coulomb 
pairing 


-0.252 
0.164 
0.065 
0.372 
0.005 

-0.020 


-0.334 
0.304 

-0.016 
0.676 

-0.161 
0.081 


total mean field 


0.334 


0.550 


correlation energy 


-0.006 


-0.051 


^ ' pn 


0.328 


0.490 


^ * pn 


0.362 


0.509 



ergy brings a large negative contribution, even for °^Er 
for which the four nuclei have very similar deformations. 
This demonstrates that the frozen approximation is not 
valid and that the use of a unique single-particle basis is 
not justified. The largest positive contribution to SVpn 
comes from the Skyrme EDF. Other terms are small for 
^^^Er, but there is a large negative contribution of the 
Coulomb term in ^^^Nd. Note also that, although the 
contribution of correlations is small for this nucleus, it 
has the right sign to bring the theoretical SVpn in good 
agreement with the data. 

Altogether, these examples indicate that, in a realistic 
model, SVpn is not determined by the interaction between 
the last two valence nuclcons, but has contributions from 
the modifications of all single-particle wave functions on 
the one hand, and from all terms in the energy functional 
on the other. 



3. Deformed rare-earth nuclei 

Let us now consider two nuclei in a deformed region 
of the nuclear chart. The four nuclei entering the calcu- 
lation of SVpn for ^®^Er have similar deformations. By 
contrast, they differ significantly for ^^^Nd, as can be 
seen in Fig. \8\ In both cases, the spherical mean-field 
results arc of no interest and will not be discussed here. 
The contribution of correlations brought by configuration 
mixing and symmetry restoration is also small and the 
analysis of SVpn can be performed for the deformed calcu- 
lations only. The value of SVpn results from contributions 
coming from all terms in the functional. The kinetic en- 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have analyzed in details the relevance 
of a difference between the binding energy of four nuclei, 
called SVpn, as a measure of the proton-neutron inter- 
action between valence particles and as an indicator for 
structural changes in nuclei. 

We have first investigated whether one can derive a 
relation between SVpn and a proton-neutron matrix el- 
ement in simple models, where analytic formula can be 
derived. Even in the oversimplified case where the four 
nuclei entering SVpn can be described by HE wave tunc- 
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tions generated by the same mean field and a two-body 
interaction, one obtains only a relation between SVpn and 
a combination of two matrix elements. Any higher-order 
term in the interaction, such as a density-dependence 
or a three-body interaction, complicates the relation, as 
do self-consistency and any correlation such as pairing, 
deformation, or any configuration mixing. This formal 
analysis already indicates that it can hardly relate 6Vpn 
to any specific proton-neutron interaction in a realistic 
model. This formal analysis is confirmed by the de- 
composition of calculated 6Vpn values for the doubly- 
magic spherical ^°^Pb nucleus and the deformed ^^^Nd 
and ^^®Er, for which we find indeed that all terms in the 
energy functional contribute, not just the proton-neutron 
interaction, and that self-consistency and correlations be- 
yond the mean field play a substantial role. 

We have then shown that our beyond-mean-field 
method has all the necessary ingredients to reproduce 
the global trends of SVpn- As has been pointed out ear- 
lier [2J], the global trend of SVpn is determined by the 
symmetry and surface symmetry energy coefficients that 
can be deduced from an energy functional. With a de- 
tailed analysis of a few representative regions in the nu- 
clear chart, we have illustrated how the fine structure 
of SVpn builds up from the successive introduction of 
deformation and correlations due to symmetry restora- 
tions and configuration mixing. Both are crucial ingre- 
dients for the reproduction of data. As found earlier for 
two-nucleon separation energies [3, [s^ , within our 
model it is essential to take into account beyond-mean- 
field correlations for the description of data in the vicin- 
ity of magic numbers. We have checked that the large- 
scale calculation of even-even nuclei using a mapped five- 
dimensional microscopic Bohr-Hamiltonian based on the 
Gogny-force [111 gives qualitatively the same results as 
ours. 

Our model provides a satisfactory description of the 
data, except at the N = Z line where the Wigner ef- 
fect is absent from our results. Within the framework 
of our model, however, SVpn docs neither provide a reli- 
able measure of the proton-neutron interaction terms in 
the energy functional used, nor a reliable indicator for 
structural changes. Certainly, in some instances struc- 
tural changes such as the onset of deformation lead to 
anomalies in the SVpn values, but in many other instances 



they do not, and there is no one-to-one correspondence 
between a structural change and the resulting modifica- 
tion of SVpn- In turn, a local increase or decrease of 
SVpn can have many different origins. One of the main 
limiting factors for the use of SVpn as an indicator for 
structural change is its asymmetric definition, which re- 
sults in the fact that additional binding from a change in 
nuclear structure contributes differently to SVpn in sign 
and size depending on the direction in the nuclear chart 
in which the structure changes. In particular the charac- 
teristic pattern of SVpn having a significantly larger size 
for nuclei where only particles or only holes are added 
to a doubly-magic nucleus as compared to systems with 
particles added for one nucleon species and holes for the 
other (a pattern recently interpreted as a "mini- Wigner 
energy" |25|). is a trivial consequence of the asymmetric 
definition of SVpm and not an indicator for a qualitative 
difference in either the proton-neutron interaction or a 
difference of their structure. 

The usefulness of SVpn is compromised by being a mass 
filter of very high order that is thereby prone to unpre- 
dictable cancelations when the nuclei entering SVpn have 
different structure. Lower-order mass filters such as two- 
nucleon separation energies or Q values are usually more 
reliable indicators of structural changes than SVpn, al- 
though these may fail as well. This also means that SVpn 
does not provide a conclusive benchmark for nuclear EDF 
methods that would be superior or complementary to 
other mass filters such as two-nucleon separation ener- 
gies and Q values. 
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